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, Breached pairing solutions to the gap equation are obtained analytically in for 

I two and three quarks and for low and high temperatures. We compare the energy 

--^ ■ of these states to that of other homogeneous states under the condition of electric 

ix. 



neutrality. We found the two-flavor BP and the three flavor mixed BCS-BP phases, 
which are stable over a wide range of parameters. Both phases contain four BP 
modes in the quasiparticle spectrum. 



.J 1 Introduction 

a \ Recently there has been considerable interest in a class of possible new states of matter 
featuring coexistence of superfluid and normal components. Examples have appeared in 
several variants, under various names ("Sarma state" P|, "interior gap" j^, "breached 
pair"0 - which we adopt here, "gapless superconductor" [7]). Breached pair states are 
candidates to arise when there are interactions favoring pairing between fermions with 
fermi surfaces of different size. They are quite different from the famous LOFF phases 
PP, and need not break translation invariance. The basic ansatz for this class of states 
goes back to old work of Sarma [H]. There are two main causes for the recent upsurge 
in interest. First, new candidate applications have emerged. These include notably cold 
atom systems, where there can be great flexibility in manipulating denisities, effective 
masses, and interactions [2/, and high-density QCD, relevant to neutron star interiors, 
as will be our concern below. Second, and importantly, parameter regimes have been 
identified where the new states are likely to be stable. General heurisitic arguments, based 
on extrapolation from clear limiting cases (ultrastrong coupling, flat bands) were presented 

*This work was completed in November, 2003 
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in j3] . Quantitative comparisons between breached pairing and other homogeneous states, 
identifying favorable cases, have been presented in [S]. 

It has also been suggested j3] that inhomogeneous (phase separated superfluid and 
normal) states can compete favorably with breached pair states. This is an interesting 
question, that receives further investigation. Our preliminary conclusion is that such 
phase separation can occur, but only in rather special parameter regimes. We shall not 
discuss it further here. 

In this paper we shall compare the energy of breached pair to conventional BCS and 
normal state alternatives, in some models suggested by high density QCD. Our calcula- 
tions may have application to neutron stars, but in view of the difficulty of interpreting 
astrophysical observations and the idealized nature of the models perhaps their main in- 
terest is methodological. We illustrate in this specific context how one establishes the 
(in)stability of breached pairing, and the infiuence of different physical conditions. 



2 Breached pairing superconductor at zero temper- 
ature 

Here we summarize properties of the breached pairing superconductor concetrating on 
the stability condition of the breached pairing phase. We present calculations in a toy 
model [S] that has a superfiuid ground state of up-strange quark pairs. The relation of 
more realistic cases to this model is worked out in the Appendix. 

We consider massive s and massless u quarks both being relativistic, 

el=P-PF^ £p=P-PF (1) 

where the fermi momenta are related to the chemical potentials as = yU„ and pp = 
\J fJ^I — fnli and the fermi momentum for the s-quark is smaller than that for the u- 
quark, p^ < p^. We suppress color indices and postulate a a weak attractive interaction 
g{,il^uO'2ips){'4'i'^2'4'l) between the light and heavy species. In a basis of light particles and 
heavy holes the quadratic part of the action is 

s-n^u (2) 

where the gap parameter is defined as A* = g /V J2p{'4^ip4'l-p) bp in a breached pairing 
(BP) superconducting ground state, and momentum sum is J2p = V I d^p/{2Ti)^. Since 
the gap parameter A is defined a s a c- nu mber, this action can be diagonalized. The 
quasiparticle energies are Spp ^ ^£p + A^, where Sp = {e^ + ep)/2 = p — pp with the 
average Fermi momentum pp = {p^ +pp)/2 and the mismatch in Fermi momenta Spp = 
(^p ~ ^p)/2 = {Pf ~ Pf)/'^ > 0- We obtain the following free energy density jS] 
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where there is no pairing in the momentum range R = {\£p\ < y^p'p ~ ^^l? which is 
singly occupied by w-quarks. The gap equation is found either by minimizing the free 
energy density or as a self-consistent condition on A 

2A r d^p A r d?p A 

(4) 



g J (27r)3 + A2 Jr (27r)3 ^^2 + ^2 

Writing / = N{0) /^^ dSp, where the density of states is A^(0) = / ^5{ep) and the 
UV cutoff is a; ~ Pi?, we have 



1 /^P.+ y/^-^-A^ 



gN{0) \AJ \ A 

Introducing the BCS gap at zero Fermi momentum mismatch, Aq, we find the breached 
pairing and the BCS solutions 

A = [Ao(25p^-Ao)]'/' {6pF>A) 

A = Ao {6pF < A) (6) 

Stability of the breached pairing state depends on the physical conditions imposed on 
the system Solving these constraints leads to the dependence of the Fermi momentum 
mismatch after pairing on the gap parameter ^^^^(A), which is parametrized [5] 

^pHA) = -^f^ (7) 

with a constant depending on a condition and 6pF is the mismatch before pairing. 
Substituting 6pF 6pf{A) in the solution of the gap equation, Eq. ©, we obtain 



A 



Ao(Ao-25j9^)/(2a2-l) (8) 



For q;2 > 1/2 one has a stable breached pairing solution 

One can check stability of the BP state also directly calculating the condensation 
energy, which is obtained by integrating the gap equation, Eq. (jH), over the gap parameter 

m 

2A' f d^p A' f d^p A' \ 

where A is a solution of the gap equation, f2jv is the energy density in the normal state. 
We use the gap equation to eliminate g, 



^BP-^N = 2N{0) ^ A'dA' In ( ^ ) - 1^ ( ^ ) + 



A'\ , /A'\ , 6pFiA') + J5pFiA')^-A 



,/2' 



5pFiA) + ^5pFiAy - A2 

(10) 
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where the first term gives the BCS condensation energy, — A^(0)A^/2, which is cancelled 
by the second term arising from integrating in the BP region R. Integrating by parts, we 
find 



Qbp ~nN = -iv(o) / 

Jo 



A A'^dA' 



d5pF{^') A' 



d^' SppiA') + ./5p^(A')2- A'2 



'11^ 



Using in Eq. (jlll) the parametrization 5p^(A) = 5pp{l + Ci'^^r), which is valid for small 

A and is equivalent to Eq. 0, we find the condensation energy of the BP state in the 
leading order O^A"^ /Spjp) 

n^p -n^ = -iv(o) , ^ (12) 

that shows a stable state when > 1/2. In the third section, we calculate for different 
BP solutions under condition of the electric neutrality, checking stability of these solutions. 

Differentiating the free energy density, Eq. (jS)), with respect to the quark chemical 
potentials, fiu and fig, we obtain the corresponding quark number densities. Namely, 
'nu + IT'S = dQpp/dpp and Uu — Ug = dQpp/dSpp, which are equal to 

f d^P (, gp \ , [ d^P gp [ d^p ^ 



Integrating, we obtain 



Pf 



n^ + ns = 2N{0) £^ , - n, = 2N{0)^6ppiAy - A^ (14) 

where the density of states is A^(0) = Pp/ir'^- When A = 0, = Pp/{'iT!'^){pp + ^Spp) ~ 
Pp^/{37i^) and Ug = Pp/ i<iT!'^)iPF — 3Spp) ~ ^^'^/(Svr^), since 6pp <^ pp, and we recover 
the free gas limit. This means that the correct way to construct the BP state is to first 
fill noninteracting quark states up to the corresponding Fermi momenta and to then pair 
that produces the condensation energy, flcond- In contrast to that, in the BCS state we 
first fill free quark states up to the common Fermi momentum, pp, and then pair 
Therefore, for the BCS state, we gain in the condensation energy, — A^(0)Aq/2, but lose 
by bringing two Fermi surfaces together, N{0)6pj^ jH]. However, the condensation energy 
of the BP state is parametrically smaller, ~ A^ Eq. (jl2j) . than that of the BCS state. 

3 Breached pairing at finite temperature 
3.1 Gap equation and quark densities 

We generalize our model to the case of nonzero temperature. The partition function at 
finite T = 1/(3 is 
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where the number of u and s-quarks, and functional integral is performed 

in momentum-frequency and in flavor spaces, Hn Dp Da c^^I n(p)'^V'a,n(p)- Assuming a 
superconducting ansatz for the ground state, in a basis of light particles and heavy holes 
the quadratic part of the action takes the form 

^ = E E K i^sA-p) ) K^, ( ^t^ji^jl) ) 



K. 



al3 



iujn -el A 
A* luj^ + el 



(16) 



where the gap parameter is defined as A* = g/VY.n^p{'^liniV')'^\ni.~V'))BP- The mo- 
mentum sum is X]p = V I d^p/{2TT)^ and the frequency sum involves Un = (2n + l)7rT. 
Integrating over Grassmann variables, we have Z = detK, where the determinant is car- 
ried out over the flavor indices and in momentum-frequency space. Employing the identity 
In deti^' = Tr In K, we find 



in^ = EEin[/3^ 



UJr, 



.fep)^ + 4 + A^ 



(17) 



n p 



where 6ep = {el—el)/2 and = {el+el)/2. Since both positive and negative frequencies 
are summed over, Eq. (fT7j) can be rewritten 



= ^ EE ( In K + W^l + A^ + 6e, 

n p 



+ In 



Following ^1 we write 



;i8) 



In [{2n + lfTT^ + (3^{ujp±6epy 
where a;, 



1 e^ + {2n + l)H^ 



+ In 



1 + (2n + 1)V 



(19) 

yp — -|- A^. The sum over n may be carried out by using the summation 
formula [Hj, ^ 

y I = Vi-^l (20) 

Integrating over 6, and dropping terms independent of (3 and fep, we finally obtain 

^ P I O/ /_9 I A9 i:_ N I /_9 , A2 



InZ 



2 

+ 2 In 



i2ny 



1 + e~^^v4+^-fep) 



el + A^-6ep)+Pi^el + A 

1 + g-/3(v/4+A^+5£p)" I 



+ 2 In 



(21) 



The thermodynamic potential is defined as = —P = —TlnZ/V. Adding the 
contribution from the mean field potential, A'^/g, and the vacuum energy arising from 
the normal ordering, ^p, we obtain 



9 
Tin 



(27 

1 + e 



Tin 



1 + e 



(22) 



where £p = p — Pr, fep = ^Pf, with pp = {pp + Pp)/2 and 6pF = {pp — Pp)/2 > 0. The 
thermodynamic potential in the four- dimensional notations, Eq. (I17j) . is 



n 



BP 



Sp 

(2 



(23) 



where TX)™ ^ = T dr = 1. Variation of the thermodynamic potential with respect to 
the gap parameter, d^Bp/dA = 0, gives the gap equation in the 3-d notations 



2A _ 
9 

rewritten as 



d^p 



A 



(24) 



2A 



d^p 



A 



9 J (2vr)3 yi^:pA^ 2 



tanh 



pUel + A2 - 5pp) 



+ tanh 



pUel + A2 + 5pp) 



(25) 



and the gap equation in the 4-d notations 



d^p 



T" (27r)3(a;„,-z5pp)2 + £2 +A2 



(26) 



where / (0 = N{<d) dsp, with the density of states N{0) = J ^S{ep) and the UV 
cutoff UJ ~ pp. 

The poles of the anomalous propagator are located at po = iun = ±yi|~+"A^ — Spp, 
which define the quasiparticle energies. They match energies at T = 0. 

Differentiating the thermodynamic potential with respect to the quark chemical po- 
tentials, Uu + rig = dflsp/dpp and riu — Ug = OVLbp / d5pp, we obtain the quark number 
densities in the 3-d 



riu + ris 



d^p 
(27r)3 

d^p 
(2^ 



1 



1 



+ 1 e 



+ I 



(27) 



and in the 4-d, 



riu + Us 



2TE 

n 

2ty: 



d^p 
(2^ 
d^p 



2 {uj,,-zSppy + el + A\ 
i{uJn - iSpp) 



(27r)3(^„-z5p^)2 + £2 



(28) 



Both the 3-d representations Eq. ()24l27p and the 4-d Eq. ()26l28p representations will 
be useful later. 
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3.2 Solving the gap equation 

Small temperatures, T <^Tc. 

We now examine how the size of the gap in the energy spectrum depends on the 
mismatch of the fermi momenta and on the temperature. 

First consider the case of low temperatures, T <C where is the critical temperature 
for the BCS superconductor, and make a suitable expansion of the gap equation, Eq. (j21|). 
We have 



gN{Q) 



de 



+ A2 



1 1 



(29) 



where A^(0) = (27r^) is the density of states, e = p~Pf, and u is the ultraviolet cutoff. 

as 

(30) 



9N{0) Ji 



dx 



Introducing a dimensionless variable x = y/e"^ + A^/A, we rewrite the gap equation as 
r/^ dx [ 1 1 

q{A/T){x+Spp/A) + 1 ~ e(A/T)(x-5pj^/A) + ^ 

In the limit A/T ^ 1 and when 6pF > A, we have 

poo rjrp °° 

Ji Vx2 - 1 

/"°° V-/_-|Nn+lp-«(A/T)(x-^pf/A) 



n=l 



^_g(A/T)(..-5pWA)^_J (31) 



When 5pp < A the third integral is over the range [1, cxo), and there is no fourth integral. 
Since the second and third integrals converge, we can set their upper limit equal to oo. 
For a not very large mismatch, Spp ~ A, we get 



in m - 2 £,-1,.- CO.. (--p.) K„ (f) . ,„ 'IZI^^ I ,32) 



n=l 



when 6pF > A, while the last term on the right-hand side is absent when Spp < A. Here, 
we introduced the BCS gap at zero temperature, Aq = Ao(T = 0), and we invoke the 
Bessel function Ko{z) = e-^'/yW^ dt 

Using the asymptotic expansion of Bessel functions at large arguments and solving 
Eq. (jSgi we find 

A(r) = [Ao(T) {25pF - Ao(r))]^/2 i5pF > A) 

Ao(T) = Ao - cosh (^^^ ^j2nTAo (l - ^) e"(^»/^) {6pF < A) (33) 

where A(T) is the BP and Aq{T) is the BCS solutions of the gap equation at low tempera- 
tures and finite Fermi momenta mismatch. At 6pF = 'we recover the known dependence 
of the BCS gap on the temperature 
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Next we calculate the condensation energy by integrating the gap equation, Eq. (piJj) . 
over the gap parameter. We use the gap equation to eliminate g, 



N 



2N{0) 



+ 2E(-ir 



+1 



n=l 



cosh 



A'dA' l^ln 
n6pF{A' 



A 



T 



In 

/nA' 



A 



T 



+ In 

cosh 



6pF{A') + J6pF{A'y~A''^ 



6pp{A) + ^6pF{Ay-A^ 
' ndpF^A) 



T 



Integrating by parts, we obtain 

r-A A'^rfA' 



BP 



n 



N 



-N{0) 







d6pF{A') 



A' 



(34) 



6pFiA'y - A' 2 

o oo 



+ iV(0) / A'2rfA'-5:(-l 

•^0 n=l 



n 



dA' 



cosh 



5pi.(A') + V'^Pf(A02- A'^ 



'n5pi7'(A') 
T 



T 



sinh 



'n6pF{A') 
T 



'nA'\ d5p^(A') 



T 



dA' 



(35) 



where we used the formula Kq{z) 



a 



2 



-Ki{z). Using the parametrization (5j9f(A) = (5pi?(l + 



25pl 



), we obtain in the leading order 0{A^ / Sp^p) 



BP 



N 



-iV(0) 



A\2a' 



S6pj, 



+ Ar(0)2T2E(-l)"+i 



n=l 



5PF 



sinh 



' ndpF^ 
T 



cosh 



f n5pF\ 



nA/T 



rfi Jo 



\ T J n"^ Jo 
K()(x)x^dx 



Ki{x)x'^dx 



(36) 



To evaluate the integrals, we write Jq^^^ Ki{x)x'^dx = 2 — /^/-p Ki{x)x'^dx, where in the 
remaining integral we use the asymptotic expansion of the function Ki (or equivalently 
we use the recursion relation for the Bessel functions / xP^^Zp{x)dx = x^^^Zp+i [IH]); and 
similar for the second integral. We obtain 

2 



nA/T 



nA/T 



Ki{x)x^dx 
Ko{x)x^dx 



2 
4 



nA 

nA-^ ^ 



K2 

K3 



nA 
~T~ 

nA 



nA 



nA 



(37) 



where we used K2{z) = IKi{z)+Kq{z), x^'K^,{ax)dx = 2^'-^a-^'-^T [^^^) T (i±f^) 
|15j . Since A/T ^ 1 at low temperatures, the Bessel functions for n = 1 give the dominant 
contribution (we use the asymptotic expansion of the function K,,{z)). For the constant 
term in Eq. (jSZj), we need to calculate the sum in Eq. In the leading order T/A <C 1, 
we obtain 



^BP — ^ 



N 



-iV(0) 



A\2a' 



86pl 



+ N{0) 



2T^ (-Li2(-e- 



(5pf/T)n 



- Li2(-e 
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-(A/T) 



-N{Q) 



6pF 



-(A/T) 



(3^ 



where the polylogarithmic functions are given by Li„(2;) = J^kLi^''/^"'^ l^^^ni^) = 
^Li„_i(z), and A is the BP solution of the gap equation given by Eq. (jHI) with Aq defined 
in Eq. When 6pF = 0, terms in the first square bracket of Eq. (jHH|) describe the 

temperature dependent part of the condensation energy in the BCS state |13j 



nBcs-nN = -Nio)^ + Nio) 



15 T 



-(Ao/T) 



(39) 



where the term ~ is the negative of the principal term in the expansion of the free 
energy of the normal state in powers of T; we used Y^kLi{~^)'^~^^ / k'^ = 7r^/12. In Eq. 
terms dependent on the gap parameter contribute to the Qbp, while both cosh and sinh 
terms give comparable contributions in the BP state when T A 6pF- This means 
that the BP heat capacity differes from that of the BCS state and is modified by imposing 
different physical conditions. 

Large temperatures, T T^. To determine the behavior of the gap for the temperatures 
near the critical temperature Tc, it is most convenient to start from the 4-d relation 
Eq. ()26|) . Near Tc the size of the gap is small, and hence in Eq. (j^UI) . we can carry out an 
expansion in powers of A^/T^ <C 1 



de 



[UJr, 



A' 



A' 



2^3 



+ ... 



(40) 



Interchanging the order of summation over the frequencies and integration over e in the 
(convergent) terms of the right-hand side, we obtain 



gN{0) 



de 



1 - 



e 

^ 2 oo 



E 



1 3 A^ °= 



(41) 



E 



i^T)' fo {{2n + 1) - ^'^y 4 (vrT)^ fo {{2n + 1) - z^) 



To evaluate the integral, we split it into two pieces, Jq dx = Jq dx + dx, where x = e/T, 
and we calculate each piece by doing necessary approximations in the corresponding 
region. Expressing the series appearing in Eq. (j4H) in terms of the Riemann zeta function. 



i.e. writing E^=o j;^- 



—C[z), we find in the leading order 



-2Ei{-V 



1 — cosh 



'5p£ 
T 



7C(3) A2 



(42) 
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where the exponential-integral function is _ Ei[—1) —0.219 

and C(3) ~ 1.202, and is the critical temperature for the BCS gap at zero mismatch. 
Thus, in the leading order, we find that the size of the gap near Tc is 



A(T) = Ao(T) 




3.06Te 



(43) 



for both the BP {5pF > A) and the BCS {6pF < A) pairing. At 6pF = we reproduce 
the known temperature dependence of the BCS gap near Tc The behavior of the 

gap as a function of the Fermi momenta mismatch for low Eq. ()33|) and high Eq. ()43p 
temperatures is shown in Fig.(l) and Fig. (2), correspondingly. To parametrize A(T) we 
used the BCS parameters (at 6pp = 0), the gap at zero temperature Aq and the critical 
temperature Tc, which are related as Tc/Aq = e'^-^/vr ^ 0.567 where 'Je ~ 0.577 is the 
Euler's constant. 

Integrating the gap equation Eq. PU)) over the gap parameter, we find the condensation 
energy density at large temperatures 



N 



+ 



2Ar(0)^^ A'dA' (2{-Ei 
7C(3) A'2- An 



-1)) 



cosh 



T 



— cosh 



'SpFiAY 



8 (7rT)2 
Integrating by parts, we obtain 



(44) 



BP 



n 



N 



-iV(0) 



7C(3)A^ 
16(7rT)2 



A^(0) / A'^dA 
Jo 



,2(-E^(-l)) 
T 



sinh 



'Spf{A')\ d5pF{A') 



T 



dA' 



(45) 



Parametrizing the Fermi momentum mismatch at large temperatures as 6pf{A) = 5pi?(l+ 
/32^), we find from Eq. (gH) 

?2/ 



N 



A^ 



7C(3) , l3\-Eti-l))6pF ^.^-^fdpF 



167r2 2 T 



7C(3) 



T 



T 



Svr^ SpF . ( SpF 



7C(3) 



sinh 



T 



+ ( 1 - - ) A{-Ez{-1)) { 1 - cosh 



'5pF_^ 



(46) 



where we used the gap equation solution, Eq. (ji^ . and kept the leading terms O^Spj^/T^). 
When 6pF = we reproduce the BCS free energy density, which is the leading term in 
Eq. (gni) since SpF <^ T^. 



10 



Next, we find the difference in the quark number densities in the breached pairing 
state at finite temperature. For low temperatures, A/T ^ 1, it is convenient to use the 
3-d formula, Eq. (|77|). 



2N{0) A 



^/^ xdx 



1 Va;^ - 1 



1 



1 



q{A/T){x-5pp/A) _^ I q{A/T){x+Spp/A) _^ I 



(47) 



where x = yiM-A^/A, and iV(0) is the density of states. As with the gap equation, 
expanding the exponentials in the integrand, we find in the leading order 



2iV(0) 



5pl - A2 



sinh 



T 



27rTAo 1 + 



ST 
^0 



-(Ao/T) 



(45 



when 5pF > A. For n„ — Ug we used integral representation of the appropriate Bessel 
function, Ki{z) = —j^Kq{z) = te~^^ /y/t"^ — 1 dt 113], and its assymptotic behavior at 
large z (see formula above). 

For temperatures near critical, we use the 4-d relation Eq. (j28|l . As with the gap 
equation, expanding in powers of A^/T^ and performing summation over n in the leading 
A-independent term and interchanging the order of summation over the frequencies and 
integration over e in the convergent terms, we obtain 



2A^(0) 



de 



1 



1 



A^ 



E 



+ ... 



i^T) „=o ((2n+l)-*^; 
Using the Riemann zeta function, we obtain in the leading order 



(49) 



2A^(0) 
2iV(0) 



rll^ + ^sinh^^^^^ 



2 T 



T 



0.50 + 0.74 T sinh 



T 



3C(2) A^ 
4 (ttT) 
A2' 

- 0.39 — 
T 



(50) 



where ({2) = 7r^/6. Sum of the quark number densities at nonzero temperature coincide 
with that at zero termperature, Eq. (HH), 



Pf 



(51) 



since the temperature dependent term in Eq. ()27|1 is an odd function of e. Particle number 
density of a free gas gets the temperature correction, Eq. for example 



Pf 
37r2 



{pf + tt'T' 



(52) 



In the following section, we use the BP solutions of the gap equation at zero, Eq. 
and nonzero temperatures, Eqs. (j33l43j) . as well as the quark number densities at T = 0, 
Eq. (IHj), and at T 7^ 0, Eqs. ()48l50l51l52p . to find the breached pairing states of the quark 
matter which are electrically neutral. 
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4 Breached pairing in electrically neutral quark mat- 
ter 



In this section we solve the gap equation(s) imposing condition of the electric neutrality. 
In addition, when we require the neutrality with respect to the color charges the energy 
minimum shifts slightly jTj. We therefore put the chemical potentials associated with 
color charges of the f/(l)3 x U{l)s subgroup to zero, fi^ = fig = 0. This approximation 
does not change our results qualitatively. 

In the color and flavor antisymmetric channel, the three flavor gap parameter is 
parametrized by three gaps, A^^, A^^, A^^, (Appendix A). One minimizes the thermody- 
namic potential fl, 

with respect to the gap parameters Aj and the electric chemical potential fie associated 
with the electric U{1) charge. This is equivalent to looking for the minimum of Q in the 
(Aj, /ie) plane only along the electric neutrality line |7j. We however do not solve directly 
minimization problem for Q. Numerical minimization was pursued for two flavors in [7j 
and for three flavors in . Here, we solve analytically the gap equations and the neutrality 
conditions and find the parameter space where they intersect. 

Expressed through the quark number densities, electric neutrality requires — 1^^^ — 
|n<j — He = 0, or written explicitly for each pair, 

1 f 2 (1) _ 1^(2) 2 _ 1 (3) _ 1^(2) _ 1 (3) 

2V3" S'^ S" 3' 3 3" 

+ f 2^(2) _ i^l) + 2 (3) _ 1 (1) _ 1 (3) _ 1 (2)\ _ (54) 

V3" 3'^ 3" 3" 3^^ 3"/ ^ ' 

where the upper index 1, 2, 3 denotes color red, green, blue, correspondingly, the first 
brackets contains quarks from the 3x3 block and the second bracket from the 2x2 
blocks (Appendix A). The quark densities are defined by the quark Fermi momenta, 
given hj pp = fi — (2/3)/ie, Pp = fi+ (l/3)/ie, = + (l/3)/ie — 't^l/'^l^ where /i is the 
baryon chemical potential and /Xe is the electric chemical potential. For convenience we 
introduce the average and the mismatch in Fermi momenta for each pair, for example, for 
Vud = {Vf + Pf)/2 and 5pud = {Pf ~ Pf)/'^ where 6p is chosen to be positive. Explicitly, 
we have 

Pud = At - TT (^Pud = — 

2 

/ie ml ml /ie 

f^e ml ^ ml , ^, 

where the ordering of the Fermi momenta in the non-interacting electrically neutral quark 
medium is PTl p'^p < p% < Pp- In what follows we specify the quark number densities 
and solve the gap equations under condition of the electric neutrality for two and three 
flavors. 
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4.1 Two and three flavor stable quark matter at zero tempera- 
ture 

Three flavor mixed BCS-BP state. At asymptotically high densities three flavor quark 
matter is in the CFL phase. Going to lower densities, the effective strange quark mass 
increases making impossible to maintain the BCS pairing for all nine quarks (the CFL 
phase). The CFL phase becomes unstable when flrst the down-strange quarks break the 
BCS pairing at m^/4yU > Aq, Eq. ()6I55|) . where Aq is the CFL gap. 

We flnd that, due to the small Fermi momentum of the strange quark, pp<^p'^< 
Pp, up blue-strange red and down blue-strange green quarks participate in the breached 
pairing while up green-down red quarks pair according to the BCS mechanism. The 
remaining up red, down green and strange blue quarks form the BCS condensate. Out 
of total nine quasiparticle excitations, there are flve BCS and four BP modes, leading to 
seven gapped and two gapless excitations (Appendix A). The system of gap equations 
for the BCS A^^ and the BP A^^ ~ A^g gap parameters decouples in the 3 channel and 
reduces to the two two-flavor gap equations (Appendix A). We therefore use the quark 
number densities and the stability criterion obtained for the two-flavor case in the previous 
sections. 

When quarks participate in the BCS pairing, their Fermi momenta and corresponding 
densities are modifled by the strong interactions in a way to equalize the densities of two 
species which form the Cooper pairs 0. In our case, n^^ = ri'P = n^J^\ meaning that 
up red, down green and strange blue quarks form an electrically neutral quark matter, 
|n^^^ — ^n^d^ + — ^nf^ — ^n^^^ — ^nf^ = 0. Therefore the electric neutrality condition, 
Eq. (|^. including quarks of other colors is given by 

- ^rid) + C^riu - - \nd. - \n)\ - rig = (56) 

J <J /BCS V<J -J <J O J BP 

where Ue = /ig/(37r^). Using Eq. (fT^ for the quark number densities, we obtain 



3 6 6 4/i 



where A = A^^- Since down and strange quarks carry the same electric charge, the 
breached region (with no pairing) from the (ds) condensate does not contribute to the 
electric neutrality. The breached region from the (us) condensate contribute to the excess 
of the positive electric charge, that already exists in the non-interacting quark matter 
where strange quarks are less abandon than up and down quarks at 7^ 0. Therefore a 
flnite density of electrons in needed, /ig 7^ 0, to satisfy the electric neutrality. 

Since Eq. (j57jl includes the gap parameter Aus, we solve the neutrality condition with 
respect to Spus = ("^s/4yU — fie/2) before (A^^ = 0) and after (A^^ 7^ 0) pairing. Using 
the parametrization 6pus{Aus) = Spusi^ + d"^^^^), we obtain in the leading order 

" 45 + 271i?2 + 24x + 231a;2 - 2i?(45 + 239a;) ^ ' 
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where x = is the numerical solution of the neutrality condition Eq. (j37|) before 

2 

pairing at some fixed strange quark mass and R = According to the stability criteria, 
Eqs. dHl??), the BP solution is stable when > 1/2. Fixing the baryon chemical potential 
H = 400 MeV and increasing the strange quark mass rris = 150 — 327 MeV, we obtain 
IPf ~ Pf\ — 14 — 56 MeV which grows and the coefficient in the stability criteria for the 
BP Aus solution = 0.59 — 0.50 which drops. Satisfying the condition of breaking the 
CFL state, Aq < ml/Afi with Aq = 20 MeV, we have m, > 179 MeV. This means that 
in this range of parameters, 179 < rris < 327 MeV and under the neutrality condition 
Eq. (jFTj) the mixed BCS-BP phase for three flavors is stable, while for larger masses, 
nis > 327 MeV it becomes unstable. 

The similar phase, containing seven gapped and two gappless excitations, was obtained 
numerically by Alford, Kouvaris, Rajagopal [F, and was called by the authors the gapless 
CFL. 

When the ordering of the Fermi momenta in the strongly interacting neutral quark 
matter follows the pattern Pp Pp < Pp, there is a hirarchi of scales between the gap 
parameters A^^ < A^^ < A^^. We might find phases containing six (when down and 
strange quarks do not pair) and eight BP modes. The phase with the maximum number 
(eight) of the BP modes requires the comparable Fermi momenta mismatches between 
different pairs and is the BP analog of the CFL phase. It seems, however, that these 
phases are not stable under condition of the electric neutality at zero temperature. 

Two flavor BP state. At large strange quark mass, only up and down quarks (up red 
with down green, and up green with down red) participate in the breached pairing, while 
strange quarks of all three colors and blue up and down quarks form the free gas. The 
electric neutrality condition, Eq. (j54j) . is written 

/2 1 \ 2 1 1 
2 -riu - -nd + -Uu - -rid - -^n^ - = (59) 
\o o / BP 6 6 6 



Using the quark number densities, we obtain 



2(/i-- 
V/- 6 



where A = A„rf. We solve the neutrality condition, Eq. (jHIIl), with respect to 5pud = l^e/2 
before (A^^ = 0) and after (A^^ ^ 0) pairing. Using the parametrization 5pud{,^ud) = 
Spudi^ + 2&v^ )' obtain in the leading order 



ud 



2 9 6x -|- X , . 

" ~ 18 - 18x + 48a;2 + 12i?2 - 4i?(3 + 2x) ^ ' 

and for the infinitely large when the free gas of the strange quarks does not participate 
in the neutrality balance, 

2 27 - 18a; + 3x'^ 

0? = 62) 

45 - 66a; + 140x2 ^ ^ 
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Increasing the strange quark mass = 200 — 400 MeV, we obtain the electric chemical 
potential /^e = 23 — 70MeV which grows and the coefficient for the stability criterion 
is = 0.53 — 0.59 which also grows. Satisfying Aq < ml/ifi with Aq = 50MeV", we 
have rus > 283 MeV. For two flavors we need higher coupling (Aq) than for three flavor 
case in order to satisfy the electric neutrality condition, i.e. the electric neutrality point 
/ie(A = 0) should be to the left from the gap equation point /ie(A = 0) = Aq (see Fig.(l)). 
For rris oo, we obtain /ig = 86 MeV and = 0.64. Therefore the BP two-flavor phase 
is stable for rris > 283 MeV, and there is no restriction on from above. 

At rris oo, this phase was obtained numerically by Shovkovy, Huang at T = [Tj 
and T 7^ [TT], [T2] and was called by the authors the gappless 2SC. 

4.2 Stability of the quark matter at nonzero temperature 

Here we find the solution of the gap equation under the condition of the electric neutrality 
at nonzero temperature graphically. We use the gap equation solutions, Eqs. ()33l43p . and 
the quark number densities, Eqs. ()48I50I51I52|) . at nonzero temperature. 



A (MeV) A (MeV) 




[ml/2fi - fi,] (MeV) fi^ (MeV) 



Figure 1: Gap as a function of the Fermi momenta mismatch, 26pF, at small temperatures. 
The BCS and BP are the two branches of the gap equation solution, and N denotes the 
electric neutrality line. Left: Three flavor mixed BCS-BP state, T = 2 MeV, Aq = 
20 MeV, rris = 2m MeV. Right: Two flavor BP state, T = 5 MeV, Aq = 50 MeV, 
rris = 300 MeV. In both cases fx = 400 MeV. 

Three flavor mixed BCS-BP state at T ^ 0. As at zero temperature, we analyze 
the BP gaps A„s ~ A^s without reference to the BCS gaps formed by the other quarks 
(Appendix A). We specify the electric neutrality for the BCS-BP state, Eq. at small 
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temperatures 



1 



(^!!^ _ ^ 



Zip m: \ 



A2 



(63) 



+ 3 sinh 



I (ml jit 



T \4fi 2 
and large temperatures 
1 




27rTAn e-(^°/^) 



2 
3 



if _ /fe _ ^ 

3 V 6 4/i 



\2T V4/i 2 y e \4fi 2 




2 / Ug m^^ ^ 

/i + — ^ 

3 V 3 4^ 



^ - /ie (/ie + TT^T^) = 



(64) 



We solve the neutrality condition with respect to A as a function of the Fermi momenta 
mismatch. In Fig.(l), the neutrality line first intersects the gap equation solution and 
then the border between the BP and the BCS states, given hj A = 6pF- Quark matter is 
positively charged to the left from the neutrality line, and it is negatively charged to the 
right. In order to get an intersection between the gap equation solution and the neutrality 
line, the Fermi momenta mismatch in the neutrality condition for a free gas should satisfy 
Spf{A = 0) < Aq where Aq is the BCS gap. We therefore have stronger coupling for the 
two flavors (Aq = 50 MeV) than for the three flavors (Aq = 20MeV). This agrees with 
the renormalization group scaling of the coupling constant as we go from high to lower 
densities. We obtain the electrically neutral solution in the range of strange quark masses 
which agrees with our stability analyses at T = 0. 



Two flavor BP state at T ^ 0. As at zero temperature, we consider the BP gap 
Aud- The electric neutrality for the two flavor BP state, Eq. ()59p . is written at small 
temperatures 



3^^^ 3 



^ ( ^ _ /f£ ) _ 3 / ( ^ ) _ A2 - 3 sinh ( ^ ) J2vrTAo e^^^"/^) 



/ie 



/ie 



(/i 



TT 



3^^ 3 



2T 



-(/i+f-i;^)((/i + t-i;;f + .^T^]-/ie(/i^ + .^T 
and at large temperatures, 

2 



2 /i- 



/ie 

6 



/i 



/ie 



3T 



H — sinh 
4T e 



/^e 

2T 



+ 



37r A2 
8"Y" 



(65) 
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Figure 2: Gap as a function of the Fermi momenta mismatch, 26pp, at large temperatures. 
Left: Three flavor mixed BCS-BP state, T = 9MeV, = llMeV, m, = 260 MeV. 
Right: Two flavor BP state, T = 24 Mel/, = 28 MeV, m, = 300 MeV. In both cases 
= 400 MeV. 



+ 1(1^ - ^) ((/i - + rr'T^) - ^(/^ + y) ((/^ + y)' + ^'^') 

+ f - if ) ((. + f - if )^ + .^T^) - (.^ + .^T^) ^ (66) 

In Fig. (2) we used the same set of parameters as in Fig.(l) (including Tc = 0.567Ao), but 
only increased the temperature. At temperatures near the critical one, the BCS and the 
BP solutions form one curve. We find that it is much simpler to satisfy the neutrality 
condition at high temperatures, that leads to a larger parameter space where the neutral 
BP phase is possible. High temperatures also open an opportunity for neutral phases 
contaning more than four BP modes. 



5 Conclusions 

We analyzed the breached pairing superconductivity at zero and finite temperatures. As 
in the previous work [21, we found that the additional constraints and physical conditions 
are crucial for the stability of the BP phase. We found analytical expressions for the 
stability criteria, showing the parameter space where the BP phase is stable under a 
general constraint. 

Imposing the condition of the electric neutrality, we found the two-flavor BP and the 
three flavor mixed BCS-BP phases, which are stable over a wide range of parameters. 
Both phases contain four BP modes in their quasiparticle spectrum. Fixing the chemical 
potential at /i = 400 MeV and increasing the effective strange quark mass, we found that 
the CFL breaks at ~ 179 MeV, followed by BP phases including the mixed BCS-BP 
phase for 179 < rris < 327 MeV and the two flavor BP phase for 283 < m^. Phases 
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containing more than four BP modes which preclude the two-flavor BP phase might be 
also possible. In the region of nis where different BP phases overlap, a phase which is 
more energetically favorable wins. 

At nonzero temperature we found solutions of the gap equation, which are consistent 
with the numerically found solutions in jHj. At low temperatures, there are two distinct 
branches, the BCS and the BP solutions, while there is only one curve for both solutions 
at high temperatures. It is much simpler to satisfy the neutrality condition at high tem- 
peratures, leading to a larger parameter space where the neutral BP phases are possible. 
High temperatures also open an opportunity for neutral phases contaning more than four 
BP modes. 
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Appendix. Gap equations for two and three flavors 

In this appendix, we obtain the gap equations for diquark condensates containing two 
and three quark flavors. We use a toy model where the full interaction between quarks 
is replaced by a four-fermion interaction with the quantum numbers of a single-gluon 
exchange, = G / (i^x(^/'(x)°A;^^7''?/'(x)f )(V^(t/)jA^57^V(l/)^), where a,P, etc. are color 
indices and i, j are the flavor indices, A"^ are the color SU{3) generators satisfying A^^^A;:^^ = 

|(3^a5^7/3 — SaisS^s)- We allow condensation in the channel A^^ = G{ilJiC'f5ipj) with the 
simplest color-flavor structure, suggested in [THI, that interpolates between the color- flavor 
locking phase at Ms = and the 2SC at large Ms, 



/ b + e b c 
b b + e c 
c c d 



/ 



/ 



/ 



(67) 



where the basis vectors are (a, i) = (1, 1), (2, 2), (3, 3), (1, 2), (2, 1), (1, 3), (3, 1), (2, 3), (3, 2) 
with colors (1, 2, 3) correspond to red, green, blue and flavors (1, 2, 3) correspond to up, 
down, strange jlHI- Detailed properties of the condensate ansatz are discussed in [TH] , 
in particular this condensate locks color and flavor. In the mean field, the four-fermion 
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interaction Li^t leads to an effective action 



Pol-E Q W ^ 
Q poI + E [^^ 



(68) 



where is a 9 component color- ffavor spinor; / 
E^f = S'^^ieJiiSji + edSi2Sj2 + e^^is^is) and Q 
condensate {l/G)Qf^Af^ , equal to 



6°'^6ij, E and Q are color- ffavor matrices, 
Lint generates also the 



a-y 



lObe + e^ + d^- 2(c^ - 6c/ + f 



(69) 



Diagonalizing the effective action matrix, Eq. ()68|) . we obtain the free energy as a sum 
of the quasiparticle energies and the condensate term. Variations of the free energy with 
respect to the gap parameters give the gap equations. Matrix is block-diagonal 
in the color-flavor space, that permits to diagonalize corresponding parts of the effec- 
tive action separately. Unitary transforming part of the effective action for the indices 



:i,l),(2,2),(3,3). 



term E^^ 



we obtain the 3x3 part of the condensate Q'^^ and the kinetic 






6 + 5e -^2(0-3/) 
-y2(c-3/) 2d 





(70) 



correspondigly, where e = \{eu + £d) and 5e = |(£m — £d)- The 2x2 parts of the are 
off-diagonal containing matrix elements (36 — e) for (1, 2), (2, 1) indices, and (3c — /) for 
(1, 3), (3, 1) and (2, 3), (3, 2) indices. 

Mixed BCS-BP three flavor state. When \eu — £d\ ^ — ^sl ~ \£d — £s\, i-e. Se ~ 0, 
system of the eigenvalue equations for the 3 x 3-block, Eq. (jHEI), decouples, producing two 
quasiparticle energies with the gap (36 — e) and four quasiparticle energies satisfying 



e^ + el + {b + 5e f + 4(c - 3/)^ + 4d^ 



(71) 



+ [e^ + (6 + 5e)2) [e^ + Ad^) + 4(c - 3/)^ (ee, - 2d{b + 5e)) + 4(c - 3/)^ = 

where the eigenvalue energy is A + po- AH quasiparticle energies with their degeneracies 
are 



Po ± V4 + (36 - e)2 (1) 

Po ± ^^sl + e^, + ib + 5e)2 + 4(c - 3/)^ + Ad^ ± ^/D (2) 

po ± ^el + {3b-ey (2) 

Po + U^u- Ss) ± ^iU^u + es)y + (3c - /)2 (2) 

Po -Ueu- Ss) ± Jaieu + es)y + {3c-f)^ (2) 



(72) 



where 
D 



+ + (6 + 5e)2 + 4(c - 3/)2 + 4rf'] 

£2 + (6 + 5e)2)(52 + u^) + 4(c - 3/)2(5,£, - 2d(6 + 5e)) + 4(c 



(73) 



3/)^ 
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We define S± = el + el + {h + 5e)2 + 4(c - 3/)2 + A<P ± /D. There are (9) x 2 
eigenvalues. In the CFL hmit (6 = c, e = /, = 6 + e, = e^), these eigenvalues reduce 
to (8) modes with the gap Ag = 3fo — e and (1) mode with the gap Ai = 8e, manifesting 
the 5'[/(3)y symmetry where V = color + flavor. Strange quark mass breaks the SU{3)v 
and 8 modes split into isomultiplets, 8 = 3 + 2 + 2 + 1. In the 3 channel these quasiparticle 
energies were obtained in ^H]. Varying the free energy, given by the sum of quasiparticle 
energies Eq. fl7^ and the condensate Eq. with respect to the 5 gap parameters, 

b, e, c, /, d: = etc., we obtain 5 gap equations. Since S± is a function of (6 + 5e) and 
(c - 3/), 5f = f and 3f = -f . Combining, 5^ - ^ = and 3^^ + ^ = 0, 
we obtain in the 3 channel (e = —b, f = — c, d = 0), 



d^p b _ ^ ^ 

\2 



(M' ^Bl + (46)= 

describing the BCS gap A^^ = 46 and the breached pairing gap A„s ~ A^^ = 4c, where 



1 f d^p c , , 

h (27r)3^(l(^^ + ^^))2 + (4c)2 



Q = { l^i^u + es)\ > ^{\{Su — £s)Y — (4c)2 } is the BP momentum integration area. 
Eq. (fTIj) contains the two-fiavor gap equations, Eq. (jl)), with G = g/A. 

BP two flavor state. When fe 7^ and fe ^ the gap parameters c = / = li = 0, 

and the eigenvalue equation for the 3x3 block reads 

- 2A2(^£2^^((36-e)2 + (6 + 5e)2) + fe2^ + (^£2^i((36-e)2 + (6 + 5e)2)-fe2 



+ 



6e^ ((36 - e) + (6 + he)f - ^ ((36 - ef - (6 + Se)'' ' 



(75) 



Solving this equation we omit terms in the square bracket, since they do not contribute 
to the gap equation in the 3 channel (e = —6). The quasiparticle energies are 



Po + fe± V£2 + i((36_e)2 + (6 + 5e)2) (1) 
7^6^r7WT76TR^ m 

(76) 



Vo-5e± ^+ |((36 - e )2 + (6 + Se)^) (1 



Po + ± + (36 - e)2 (1) 



which are (4) x 2 modes. Condensate is ^2 (6^ + 106e + e^). Combining, 3^ + ^ = 0, 
and taking the 3 channel, we obtain 

^^ + 2/^^^ = (77) 
G Jq (27r)3 ^^2 + (45)2 



describing the breached pairing gap A^^ = 46, where Q = {\e\ > ^JSe'^ — (46) ^ } is the 
BP momentum integration area. 
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